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Interaction between collective monopole oscillations of a 
trapped Bose-Einstein condensate and thermal excitations is 
investigated by means of perturbation theory. We assume 
spherical symmetry to calculate the matrix elements by solv- 
ing the linearized Gross-Pitaevskii equations. We use them to 
study the resonances of the condensate induced by tempera- 
ture when an external perturbation of the trapping frequency 
is applied and to calculate the Landau damping of the oscil- 
lations. 
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I. INTRODUCTION 

Since the discovery of Bose-Einstein condensation in 
magnetically trapped Bose gases, the study of the low- 
energy collective excitations has attracted a big interest 
both from experimental and theoretical point of view. 
Mean field theory has proven to be a good framework 
to study static, dynamic and thermodynamic properties 
of these trapped gases. In particular, it provides pre- 
dictions of the frequencies of collective excitations that 
very well agree with the observed ones. Recently [E10] 
the energy shifts and damping rates of these low-lying 
collective excitations have been measured as a function 
of temperature. However, these phenomena have not yet 
been completely understood theoretically. 

In this paper we study the influence of thermal exci- 
tations on collective oscillations of the condensate in the 
collisionless regime. Previous papers on this subject has 
been devoted mainly to calculation of the Landau damp- 
ing by means of perturbation theory. In Refs. 0-0] only 
the uniform system has been considered, whereas in Refs. 
JJJU Landau damping in trapped Bose gases has also 
been studied but using the semiclassical approximation 
for thermal excitations and the hydrodynamic approxi- 
mation for collective oscillations. An important point of 
Refs. 0,|| is that the authors discuss the possible chaotic 
behavior of the excitations in an anisotropic trap. The 
frequency shift has also been studied for a trapped con- 
densate in the collisionless regime in Ref. M. 

In the present work we study the interaction be- 
tween collective and thermal excitations using the Gross- 
Pitaevskii equation and perturbation theory. We con- 
sider spherically symmetric traps, since in this case the 



spectrum of excitations is easily calculated, avoiding the 
use of further approximations. Even though the case 
of anisotropic traps can be significantly different in the 
final results, a detailed investigation of spherical traps 
is instructive. We explore, in particular, the properties 
of monopole oscillations by studying the temperature- 
induced resonances that occur in the condensate when 
an external perturbation of the trapping frequency is ap- 
plied and, also, the Landau damping associated with the 
interaction with thermally excited states. 

This paper is organized as follows. In section II we in- 
troduce the general equations that describe the elemen- 
tary excitations of the condensate within the Bogoliubov 
theory M. In Sec. Ill we recall the perturbation theory 
for a trapped Bose-condensed gas in order to study the 
interaction between elementary excitations. In Sec. IV 
we introduce the linear response function formalism and 
calculate the response function of the condensate when a 
small perturbation of the trapping frequency is applied. 
We derive analytic equations for the response function 
at zero temperature and treat perturbatively the contri- 
bution of the elementary excitations, which is related to 
Landau damping. In Sec. V we discuss the main results. 



II. ELEMENTARY EXCITATIONS OF AN 
ISOTROPIC TRAP 

We consider a weakly interacting Bose-condensed gas 
confined in an external potential V ex t at T = 0. The 
elementary excitations of a degenerate Bose gas are as- 
sociated with the fluctuations of the condensate. At low 
temperature they are described by the time dependent 
Gross-Pitaevskii (GP) equation for the order parameter 
00: 



m|*(r,*) = 



2m 



+ Mr)+ fl |tt(r,i)ntt(r,i), 

(1) 



where J"dr | \J/ 1 2 = Nq is the number of atoms in the 
condensate. At zero temperature it coincides with the 
total number of atoms N, except for a very small dif- 
ference SN -C N due to the quantum depletion of the 
condensate. The coupling constant g is proportional to 
the s-wave scattering length a through g = AnTi a/m. 
In the present work we will discuss the case of positive 



scattering length, as for 87 Rb atoms. The trap is in- 
cluded through 14 x t, which is chosen here in the form of 
an isotropic harmonic potential: V cxt (r) = (l/2)raw^ r 2 . 
The harmonic trap provides a typical length scale for the 
system, a ho = (h/mto^o) 1 / 2 . 

So far experimental traps have axial symmetry, with 
different radial and axial frequencies, but experiments 
with spherical traps are also feasible |Q. The choice 
here of a spherical trap has two different reasons. First, 
it greatly reduces the numerical effort and will allow us to 
study the interaction of oscillations with elementary exci- 
tations without any further approximations. Second, the 
energy spectrum of the excitations in such a trap is well 
resolved yielding to the appearance of well-separated res- 
onances. In anisotropic traps, conversely, the spectrum 
of excitations is much denser. 

The normal modes of the condensate can be found by 
linearizing equation (Ffl) , i.e., looking for solutions of the 
form 



*(r,t) 



_ g-^A 



*o(r)+w(r)e-^'+w*(r)e i 



(2) 



where p is the chemical potential and functions u and 
v are the "particle" and "hole" components characteriz- 
ing the Bogoliubov transformations. After inserting in 
Eq. (pi) and retaining terms up to first order in u and v, 
one finds three equations. The first one is the nonlinear 
equation for the order parameter of the ground state, 



(flo + S*o(r))*o(r)=/i*o(r), 



(3) 



where Ho = — (H /2m) V 2 + T4xt( r ); while u(r) and v(r) 
obey the following coupled equations [O] : 



Tluju(y) - 


= [H - 


-p + 2g^ 2 ]u{r)+g^ 2 v(r) 


(4) 


hjjv{v) - 


= [H - 


-p + 2g^ 2 }v{r)+g^ 2 u{r). 


(5) 



Numerical solutions of these equations have been found 
by different authors |l^-|lq]. In the present work, we use 
them to calculate the response function of the condensate 
under an external perturbation and the Landau damping 
of collective modes. 

When the adimensional parameter Na/a-t is large, 
the time-dependent GP equation reduces to the hydro- 
dynamic equations |l9[: 



m m v 



dp 
at 



V(K, 



V(vp) = 



gp 



-) = o, 



(6) 



(J) 



2 is the particle density and the 
(**V*-*V**)?i/(2mip). The 



where p(r,t) = | ^(r, t) 
velocity field is v(r, t) = 

static solution of equations (ph-(ffl) gives the Thomas- 
Fermi ground state density, which in the spherical sym- 
metric trap reads 



p(r) = <rV-Kxt(r)] 



(8) 



in the region where p > V cx t(r), and p — elsewhere. 
The chemical potential p is fixed by the normalization of 
the density to the number of particles Nq in the conden- 
sate. The density profile (||) has the form of an inverted 
parabola, which vanishes at the classical turning point 
R defined by the condition p = V cx t(R). For a spherical 
trap, this implies 



p 



mujZ„R 2 



■&- 



(9) 



It has been shown Jl9[ | that the hydrodynamic equa- 
tions (||) and (0) correctly reproduce the low-lying nor- 
mal modes of the trapped gas in the linear regime when 
Na/a-to is large (see however Ref. p0|). 



III. PERTURBATION THEORY 

Let us briefly recall the perturbation theory for the in- 
teraction between collective modes of a condensate and 
thermal excitations as it was developed in Ref. Q. Sup- 
pose that a certain mode of the condensate has been ex- 
cited and, therefore, it oscillates with the corresponding 
frequency f2 sc- We assume that this oscillation is clas- 
sical, i.e. the number of quanta of oscillation (n osc ) is 
very large. Then, the energy of the system associated 
with the occurrence of this classical oscillation can be 
calculated as E — hfl osc n osc with n osc 3> 1 . Due to in- 
teraction effects, the thermal bath can either absorb or 
emit quanta of this mode producing a damping of the 
collective oscillation. The energy loss can be written as 



E = -nn osc (w^-w^), 



(10) 



where W^ a > and W^' are the probabilities of absorption 
and emission of one quantum hQ osc , respectively. The 
interaction between excitations is small, so one can use 
perturbation theory to calculate the probabilities for the 
transition between a ?-th excitation and a A:-th one, avail- 
able by thermal activation 



W = TT 



Z^ 



I (k I V int | i) | 2 



(11) 



Let Ei and E^ be the corresponding energies and assume 
Ek > Ei. Since energy is conserved during the transition 
process, one has Ef. — Ei + %Ct osc . 

The interaction term in second quantization is given 
by 



V- nt = - / </r*tvT/ T vl/vl/. 



(12) 



In the framework of Bogoliubov theory, the field opera- 
tor \l/ can be written as the sum of the condensate wave 
function ^o, which is the order parameter at equilibrium, 
and its fluctuations <5*, where * = * + £* [see Eq. (S)]. 



The fluctuations can be expressed in terms of the annihi- 
lation (a) and creation (or) operators of the elementary 
excitations of the system: 



*# = £)[«i(r)aj 



v*{v)a]] 



(13) 



where the functions u and v are properlynormalized solu- 
tions of equations (H)-(@)- In the sum ( [l3[ ) one can select 
a low energy collective mode, for which we use the no- 
tation u osc , v osc , a osc , ctJsc) an( i investigate its interaction 
with higher energy single-particle excitations, for which 
we use the indices i, k as in (|ll]). These latter excitations 
are assumed to be thermally excited. Inserting expres- 
sion (H3) into Eq. (O) one rewrites the interaction term 
Mnt in terms of the annihilation and creation operators. 
Since we want to study the decay process in which a 
quantum of oscillation hQ osc is annihilated (created) and 
the z-th excitation is transformed into the k-th one (or 
viceversa), we will keep only terms linear in a osc (aj sc ) 
and in the product a k ai (akCq). And the energy conser- 
vation during the transition process will be ensured by 
the delta function 5(E k — E% — fr^osc)- This mechanism 
is known as Landau damping pl| . 

Assuming that at equilibrium the states i, k are 
thermally occupied with the usual Bose factor /j = 
[exp(EifkBT) — l] -1 , the rate of energy loss can be cal- 
culated as 101 



E 



E = -27T- J2 I A ik I 2 HEk - E t - hQ osc )(fi - f k ) 



ik 



(14) 



where 



Ak = 2.9 / drip [(utvi + v%Vi + u* k Ui)u osc 

+ {v* k Ui + v* k Vi + u* k Ui)v osc }. (15) 

Let us define the dissipation rate 7 through the follow- 
ing relation between the energy of the system E and its 
dissipation E: 



E = -2-yE. 
Using expression ( |l4| ) 7 can be calculated as 

= ^2 ~fik s(ujik - ri osc ) , 



a, 



(16) 



(17) 



ik 



where the transition frequencies ujik = {E k — Ei)/h are 
positive. The "damping strength'' 



lik 



n 2 n a 



I A ik I 2 (fi - f k ) 



(18) 



has the dimensions of a frequency. In this work we cal- 
culate the quantities 7^ by using the numerical solutions 
u and v of Eqs. (||-§|) m to the integrals (|l^). The results 
will be discussed in section V. 



IV. RESPONSE FUNCTION 

The results of the previous section can be used also 
to study the effect that an external perturbation of the 
trap has on the collective excitations of the condensate. 
Let us assume the trapping frequency in the form [who + 
<5who(£)L where &j n o ~ exp(— iu>t) is a time-dependent 
modulation. Assuming that the perturbation is small, 
one can use the response function formalism to describe 
the fluctuations of the system. Let us briefly recall the 
basic formalism p2j| . 

The behaviour of a system under an external pertur- 
bation can be described by studying the fluctuations that 
may generate the external interaction to a certain phys- 
ical quantity of the system. An external perturbation 
acting on the system is described by a new term in the 
Hamiltonian of the type 



V = -xf(t) 



(19) 



where x is the quantum operator of the physical quantity 
that may fluctuate, and f(t) is the "perturbing force". 
The mean value (x) is zero in the equilibrium state, in ab- 
sence of perturbation, and is not zero when it is present. 
For a periodic perturbation f(t) ~ exp(— iu>t), the rela- 
tion between (x) and f(uj) is 



(x) =a(u)f , 



(20) 



where a(u>) is the response function also called gener- 
alised susceptibility. 

In general a{u>) is a complex function. It can be seen 
that the imaginary part of the susceptibility determines 
the absorption of energy Q of the external force / by the 
system through the following relation: 



Q = -Im[a(u;)]\f\\ 



(21) 



and that the real and imaginary parts of a(to) satisfy the 
Kramers-Kronig relation 



Re[a(u>)] 



P 

7T JO 



^ — OJ* 



(22) 



where P means the principal value of the integral. 

The time-dependent external drive <5w n o induces oscil- 
lations of the condensate density Sp with frequency lo; 
p{r, t) = p(r, 0) + Sp. Expanding the energy due to the 
confining potential, E\ w = fV C xtpdr, with respect to 
Sujh and Sp one obtains the "mixed" term, correspond- 
ing to the Hamiltonian (J19|): 



V = mwhoSiOho / r 2 5pdv . 



(23) 



Comparing it with Eq. (M) one can identify the perturb- 
ing force and the corresponding coordinate as 



/ = -muJhoSoJho , x= I r 2 Sp(r,t)dr 



(24) 



Note that the first order term mwho^Who / r 2 p(r, 0)dr can 
be omitted because it gives an additive shift in the Hamil- 
tonian which does not contribute to the equations of mo- 
tion of the system. 

Once we have identified / and x, we can calculate the 
response function of the condensate a(u>). According to 
the definition one has 



a(w)f. 



(25) 



Let us present the response function in the form a(w) = 
<xq(w) + ai(u)), where ao(u) corresponds to the response 
function of the condensate at T = 0, i.e., calculated with- 
out elementary excitations, and ot\(oj) is the contribution 
of the excitations. At low temperatures it can be as- 
sumed that ai(u>) <C ao(ui) and then a\ can be treated 
as a perturbation. 

We proceed as follows. First, we use the hydrodynamic 
approximation to obtain the response function at T = 
0. Then, within a perturbation theory, we introduce the 
contribution of the elementary excitations at finite T to 
obtain a\(ui). 



A. Calculation of ao(w) at T = 

For a spherically symmetric breathing mode fl23| , one 
can easily prove that the hydrodynamic equations of mo- 
tion (0) and (m) admit analytic solutions of the form E4] 



p(r, t) = a (t) - a r (t)r 2 , v(r, t) = a r (t)r . 



(26) 



These equations are restricted to the region where p > 0. 
Notice that they include the ground state solution (g) in 
the Thomas-Fermi limit. This is recovered by putting 



0. 



mu>^ /(2g), and ag = p/g. Insert- 



ing Eqs. (EG) into the hydrodynamic equations, one ob- 
tains two coupled differential equations for the time de- 
pendent coefficients a r (t) and a r (t), while at any time 

a = —(15N/8tr) 2 / 5 a r ' is fixed by the normalization of 
the density to the total number of atoms. The form (EG) 
for the density and velocity distributions is equivalent to 
a scaling transformation of the order parameter. That is, 
at each time, the parabolic shape of the density is pre- 
served, while the classical radius R, where the density 
(p6|) vanishes, scales in time as pj| 



R{t) =R(0)b(t) 



2/y 



b(t) 



(27) 



where the unperturbed radius i?(0) is given by Eq. fl9). 
The relation between the scaling parameter b(t) and 



the coefficient a r (t) is a r 
into Eq. (pi) we obtain 



"ho 



/(2gb 5 ). Inserting it 



P(r,t) 



1 



1 



^W 



26 5 ^b 3 



Pi 



(28) 



And the hydrodynamic equations then yield a r = b/b 
and 



[u ho + Suj ho (t)] 2 b 



6 4 



0. 



(29) 



The second and third terms of (129) give the effect of the 
external trap and of the interatomic forces, respectively. 
From (E7h and (E3) it follows that at equilibrium b = 1 
and 6 = 0. For a small driving strength <5who, one can 
assume that the radius of the cloud is perturbed around 
its equilibrium value, so 



R(t) = R(0) + SR(t) , b(t) = 1 + 5b{t) , 



where 



Sb{t) 



5R(t) 



(30) 



(31) 



It means that 8b is the fractional amplitude of oscillations 
of the radius and, therefore, it is a measurable quantity. 
In the small amplitude limit, one can linearize Eq. (|29j) 
with respect to i5who and 8b yielding the following equa- 
tion 



8b + 5ui^ 8b = — 2who&^ho . 



The solution is 



Sb(t) 



—2u>ho 

li M — U) 



8ui 



ho ■ 



(32) 



(33) 



where Om = v5 Who corresponds to the frequency of the 
normal mode of monopole in the hydrodynamic limit p9[ . 
Keeping only the lowest order in the small perturbation 
8b, Eq. (Bn) yields 



p(r, t) = p(r, 0) + 



c 2 2 o 



8b, 



(34) 



and using the Thomas-Fermi radius at equilibrium 
it follows that the density fluctuation is given by 



Sb(t) . (35) 

and (pa), finding 



dp(r,t) = p[r, 


t) - p{r, 0) 


5/i 
9 


U(0)J 5 


We can calculate now x using (|j) £ 


X 


= C8b{t), 



where C = 16irpR(0) 5 /(35g). Then, from Eq. (| 
gets 



-2Cwho 
7)2 772" 



Sujh 



(36) 



one 



(37) 



At T = there are no thermally excited states and, 
hence, a{ui) = ao{ui). By comparing the definition ( pa ) 
with (P7l) one has 



a {uj) 



2C 



m(0^ - lu 2 ) 



(38) 



This is the response function at zero temperature without 
including any dissipation. Therefore ao(uj) is real, i.e., 
the induced oscillations at T = are undamped. 

The energy of oscillation can be calculated as twice 
the mean kinetic energy associated to the mode, E — 
J dr p(r, 0)i> 2 . For a monopole mode in an isotropic trap, 

the calculation g|,Ei § ives E = T N lA 5b \ 2 i where \ Sb \ is 
the amplitude of the oscillation of the cloud ( |3l| ) . Using 
Eqs. <M) and S at T = 0, it follows that 



E 



15 
7G 2 



fiN\a (oj)f\' 



(39) 



B. Calculation of qi(w) 

Now we want to calculate the contribution of the ther- 
mally excited states to the response function. We study 
the low temperature regime, where u\ <C «o ancl the 
energy of oscillation ([39]) can be estimated using ao in- 
stead of a(io) — ao(u>) + ai(w). The effect of a± will be 
introduced within a perturbation theory. 

We have already seen that the thermal excitations can 
either absorb or emit quanta of oscillation fito and thus 
they will dissipate energy. The contribution of the ele- 
mentary excitations to the susceptibility will be a com- 
plex function, a\{u>) = Re[a{\ + ilm[ai], whose imag- 
inary part is related to the absorption of energy Q of 
the external perturbation. However, in a stationary solu- 
tion which is the case under consideration, the absorption 
Q must be compensated by the energy dissipation ( |l6| ) 
due to the interaction with the elementary excitations. 
Therefore, 



(40) 



E = 0. 



Let us rewrite the definition of the damping rate (Iq) by 
using ([IT]) and (|l8|) with a generic oscillation frequency 

lu: 



E = -2oj ^2 jik S(uj lk - w)E . 



(41) 



ik 



Inserting Eq. 



and defining /3(uj) = ao(u>)/C 



2/ [m(f2jy[ — uj 2 )] one obtains the energy dissipation 

jB = _2c^^ 7lfc ^ fc - W )|/3H| 2 |/| 2 . (42) 

ik 

Let us recall that the energy dissipation according to 
Eqs. ( pi|) and ( f40|) can be calculated also from the imagi- 
nary part of the response function a(cu) = ao(uj) + a\(uj) . 
Since ao(w) is real, Eq. ( pl| ) becomes 



= -Im[a 1 {uu)]\tf = -E. 



(43) 



Comparing Eqs. (|42j) and ([13|) one can calculate the imag- 
inary part of ai(uj) as 



Im[ai(aj)] 



— — ^2 ^ lk S ^ lk ~ w ) \Pi u ) i 2 ■ 



7 



(44) 



ik 



And using the Kramers-Kronig relation (E2) one finds the 
real part 



Re[ ai (co)] = -J2 



oJikjik 15 ^N 



*^u 2 k -oj 2 7 



(^fc)l 2 . (45) 



Now we have all the ingredients to calculate the response 
function of a spherically symmetric trapped condensate 
when the monopole mode is excited and a small pertur- 
bation of the trapping frequency Slu^ ~ exp(— iuit) is 
applied. It can be calculated within first order pertur- 
bation as a(ui) = ao(u>) + Re[ai(uj)] + \Im[ai{uj)], by 
using Eqs. (|38|), (E^) and (H), respectively. It is worth 
stressing that the real part of the susceptibility diverges 
at £1m (resonance of the condensate at T = 0) but also 
at ojik which are the frequencies of the thermal excited 
modes that due to the interaction are coupled with the 
monopole. 

Actually, the resonances of the condensate can be 
found by measuring the fractional amplitude of oscilla- 
tions of the cloud radius 5b at different perturbing fre- 
quencies. This measurable quantity can be easily related 
to the response function a(ui) from Eqs. (Bq) and (E5ft 



m£! ho 
-a{w) (Who • 



(46) 



Note that the perturbation theory we have used is valid 
when | u\ |<c| ao |. This condition becomes very restric- 
tive at lu near u>m- However it is not difficult to improve 
the approximation in this region by taking benefit of the 
analogy between the response function and the Green 
function G. 

It is well known that the Green function obeys the 
Dyson equation [£6| which relates the perturbed quantity 
(G) and the unperturbed one (Go) through the inversed 
functions (G _1 and Gq ) in such a way that a perturba- 
tion theory for G _1 has a more wide applicability that 
for G. Analogously, we will find a relation between the 
inverse response functions, perturbed (a -1 ) and unper- 
turbed (ttp ). One has 



1 
a 



1 



1 



(ao + ai) ao(l + ai/ao) 
and formally with the same accuracy 



(47) 



a ao ao 

m 



— (n^ -«')-- 



ik 



7G 2 



0J7. 



■CJ* 



(48) 



Now the applicability of (|4S]) is restricted only by the 
condition that the second term is small with compare to 

_m_(~)2 

2C M- 



It is worth noting that according to the equation pq ) 
the poles of a(w) related to the resonances are shifted 
with compare to frequencies u>ik and are given by the 
equation: a^u'p) / a^{uj' R ) = 1. However these shifts are 
very small. 



V. RESULTS 

In order to present numerical results we choose a gas 
of 87 Rb atoms (scattering length a = 5.82 • 10~ 7 cm). 
For the spherical trap we fix the frequency Who = 2-7T187 
Hz, which is the geometric average of the axial and radial 
frequencies of Ref . |Tj , and corresponds to the oscillator 
length aho = 0.791 • 10~ 4 cm. We solve the linearized 
Gross-Pitaevskii equations Eqs. (||-|j) at zero tempera- 
ture, to obtain the ground state wave function ^o an d 
the spectrum of excited states Ei as well as the corre- 
sponding functions Uj (r) , Vi (r) . In spherically symmet- 
ric traps the eigenfunctions are labeled by i — (n,l,m), 
where n is the number of nodes in the radial solution, I 
is the orbital angular momentum and m its projection. 
The eigenfunctions are u n i m {v) = U n i(r)Yi m (9,ip), the 
energies E n i are (21 + 1) degenerate and the occupation 
of the thermally excited states is fixed by the Bose factor. 

For a fixed number of trapped atoms, N, the number of 
atoms in the condensate, A , depends on temperature T. 
At zero temperature all the atoms are in the condensate, 
except a negligible quantum depletion Eg ] . At finite tem- 
perature the condensate atoms coexists with the thermal 
bath. In the thermodynamic limit |27f] the T-dependence 
of the condensate fraction is N (T) = N[l - (T/T c ) 3 ]. 

We consider the collective excitations in the collision- 
less regime. This regime is achieved at low enough tem- 
perature. The excitation spectrum at low temperature 
can be safely calculated by neglecting the coupling be- 
tween the condensate and thermal atoms pa ]. It means 
that the excitation energies at a given T can be obtained 
within Bogoliubov theory at T = normalizing the num- 
ber of condensate atoms to Nq(T). 

We investigate the monopole mode (I = m = and 
n = 1). The functions w osc and u osc do not present angu- 
lar dependence, and from Eq. ( |lq ) it is straightforward 
to see that the matrix element Aih couples only those 
energy levels (i, k) with the same quantum numbers I 
and m. That is, the selection rules corresponding to the 
monopole- like transition are Al — and Am = 0. It is 
obvious, also, that different pairs of levels with the same 
quantum numbers n and I but different m give the same 
contribution. Therefore, only the integration of the ra- 
dial part has to be done numerically. 

Fixed Nq and at a given temperature, we calculate the 
damping strengths dlq) for the transitions uo^ coupled 
with the monopole. In Figure 1 we show the values of 
"fik (in units of the frequency of the monopole Clyi) for 
N = 50000 87 Rb atoms at k B T = fj,. The arrow points 
to the frequency of the breathing mode Dm = 2.231 Who, 



and the chemical potential is fi = 15.69 hio^o [these values 
are numerical results of the linearized Gross-Pitaevskii 
equations (||)-(||) for N = 50000 rubidium atoms]. The 
position of the bars correspond to the allowed transition 
frequencies w^ (in units of Who) whereas their height de- 
fines the numerical value of 7^ J29[ . 

One can see that there are two different types of al- 
lowed transitions u>ik- The damping strength associated 
to most of them is very small. Conversely, there are a few 
transitions which give relatively large values of 7^ . The 
latter correspond to transitions between the lowest levels 
(rik = l,7ij = 0) for different values of I (I = 2,3,4,5). 
The main reason for these "strong transitions" is that the 
temperature occupation factor for these low-lying levels is 
large. Moreover the calculation shows that the matrix el- 
ements are also enhanced compared to other transitions. 
This is due to the fact that the radial wave functions in- 
volved in the integration have either one (rik = 1) or no 
node (rii = 0), differently from the oscillating character 
of the radial wave functions associated to higher levels 

& 

The contribution of the other transitions is like a small 
"background" which is difficult to resolve in the scale of 
the Figure. A close-up view of the damping strengths 
of the transition frequencies around the monopole is dis- 
played in the inset of Fig. 1 in order to show the dense 
background. It is worth stressing that such a distinction 
between "background" and "strong" transitions depends 
on the number of condensed atoms in the system and, 
of course, on temperature. When the number of atoms 
in the condensate increases, the number of excited states 
available by thermal excitations also increases, leading to 
a denser and less resoluble background. 

In Figure 2 we present the same as in Fig. 1 but for 
A = 5000 atoms of rubidium at k B T = fi, where here 
fi = 6.25?iwho. In this case, one can see that the difference 
between the "strong" and "weak" transitions is not so 
impressive as in a bigger condensate since all damping 
strengths can be appreciated in the same scale. 

We can conclude that at large Ao we have actually 
two different phenomena. The strong transitions create 
temperature induced resonances which can be observed 
in direct experiments. The background transitions give 
rise to Landau damping of the collective oscillations (See 
subsection B). 



A. Temperature-induced resonances 

Using the transition frequencies LOik and the corre- 
sponding damping strength 7^., we have calculated the 
response function a(ui). At zero temperature, the re- 
sponse function ocq(ui) given by Eq. (J3q), gives a reso- 
nance at the monopole frequency J7m = v5who evaluated 
in the hydrodynamic regime. Due to interaction, ther- 
mal excited modes are coupled with the monopole. It 
means that when one excites the breathing mode of the 



condensate, the elementary excitations can give rise to 
other resonances at u>ik, which are the frequencies where 
Re[a.\{uj)\ diverges [see Eq. (Bq)]. We will now discuss 
the conditions for the observation of these effects in ac- 
tual experiments. In particular, we calculate the contri- 
bution of these resonances to the response function and 
estimate the associated strengths. 

Let us study the resonances at fc^T = /i for JVo = 
150000 atoms of 87 Rb. The behavior of the damping co- 
efficients -fa- is analogous to the one for 50000 condensate 
atoms (see Fig. 1) but in this case the difference between 
" strong" resonances and small background is even big- 
ger: the dense background is not more resoluble in the 
scale of the strong resonances. There are five resonances 
that stand out the others, and that we label as uir and 
7ij the corresponding damping strength (see table 1 for 
numerical values). 

For perturbing frequencies close to the monopole 10 ~ 
Hm: the monopole susceptibility , Eq. (pq), can be ap- 
proximated to 



a Q (uj) 



2C 



V 



1 



m(flu — w)(Om + w) (Qm — w) 



(49) 



where Aq = C /(m£lyi). 

Analogously, a\(u>) near each resonance u ~ lor can 
be presented in the form ai(u>) ~ A\/(uir — to). The 
ratio Ai/Aq is a measure of the relative intensity between 
temperature-induced and monopole resonance. 

Table 1 displays the numerical values of the relative in- 
tensity for each temperature-induced resonance lor with 
respect to the monopole one, for Ao = 150000 atoms in 
the condensate at fc^T = /i. The relative strenght of the 
response function (Ai/Aq) at lor depends not only on 



the damping coefficient 7^ but also on (O 



i)- 1 - It 



means that one mode ujr will be easier to excite, i.e., the 
strength of the response will be bigger, when it is close 
to the frequency of the monopole. Note also that the 
resonance strength increases with temperature through 

From table 1 one can see that the biggest resonance 
occurs at lur — 2.2576 Who which is resoluble from the 
monopole frequency 17m = 2.234 Who and has a large 
enough relative strength to be observed. It means that, 
tuning the perturbation frequency w to this value, a fluc- 
tuation of the fractional amplitude of oscillations can be 
observed. 

In Figure 3 we have plotted the frequency dependence 
of the real part of response function a(w) calculated ac- 
cording to equation ( |4q ) for No = 150000. The response 
function is given in arbitrary units, and frequency is in 
units of Who- The dashed line shows the monopole reso- 
nance at 12m, whereas the other divergences of a(w) cor- 
respond to the temperature- induced resonances at lor. 
From this figure one can see that the thermal induced 
resonances are quite distinct one from the other and 
from the monopole one. Therefore, temperature-induced 
resonances could be observed in experiments with good 



enough frequency resolution and good accuracy in the 
measurement of the radius fluctuations. 

We would like to stress that the phenomenon we have 
discussed is related to quite delicate features of interac- 
tion between elementary excitations, and therefore, its 
observation would give rich information about properties 
of Bose-Einstein condensed gases at finite temperature. 



B. Landau damping of collective modes 

From Fig. 1 one can see that the weak background 
transitions u>ik have, generally speaking, very small fre- 
quency separation. To estimate this distance quantita- 
tively let us renominate the resonances by an index i in 
the order of increasing value of w. Then, one can define 
the average distance between resonances Aw according 
to: 



_ = £. 7 ,(w i+1 



Wj) 



Ei7i 



(50) 



In a small interval around the collective oscillation 
0.82 J7m < w^ < 1.18 Om, we sum up all the transi- 
tion frequencies allowed by the monopole selection rules 
and find the following values for the average distance be- 
tween two consecutive transition frequencies: Aw/who — 
0.0006,0.001 and 0.006 for N = 150000, 50000 and 
5000, respectively. It is hopeless, of course, to try to 
resolve these resonances. Actually, there are reasons to 
believe that these resonances are smoothed and over- 
lapped. First of all, because a real trap cannot be ex- 
actly isotropic. It means that levels with different m have 
slightly different energies, only levels with m = ± | m \ 
are exactly degenerated. Therefore, each energy level 
with a given I will be splitted on I + 1 closer sublevels 
making more dense the energy spectrum. Furthermore, 
all excitations at finite temperature have associated a fi- 
nite life time. Excitations with E ~ /x which are the 
ones that mainly contribute in the "background" tran- 
sitions, have the shortest life time. This can be ac- 
counted for phenomenologically by assuming that these 
levels have a finite lorentzian width A. That is, in- 
stead of delta functions in the equation for the damp- 
ing rate (p/7|), we will consider a Lorentzian distribution 
centered at u>ik with a fixed width A: /i(wifc,A) = 
A/(27r?i)[(a;iA: — ^ osc ) 2 + A 2 /4]. In this case the damp- 
ing rate becomes a smooth function of Q osc and its value 
when fi osc = Om defines the Landau damping of the 
monopole oscillations. At conditions 



Aw < A < w h( 



(51) 



the damping rate will have only a weak dependence on 
the exact value of A. In Figure 4 we plot the the di- 
mensionless damping rate 7/JIm as a function of the 
lorentzian width A (in units of Who) for Ao = 50000 at 
different temperatures. The summation in (flTn has been 



done over all resonances excluding of course the "strong 
resonances" presented in Figures 1, 2 and 3. One can 
see that the A-dependence is weak indeed in the interval 
A/o^ho = 0.05 -j- 0.2 and 7 can be reliable extrapolated 
from this interval to the value A = 0. We take as Landau 
damping this extrapolated value of 7. One can estimate 
the accuracy of this extrapolation procedure to be of the 
order of 10 % according to the change of 7 over this inter- 
val. In Figure 5 we plot the damping rate versus fcsT//z 
for No — 150000 and 50000 atoms in the condensate. 
As expected, Landau damping increases with tempera- 
ture since the number of excitations available at thermal 
equilibrium is larger when T increases. One can distin- 
guish two different regimes in Fig. 5 one at very low T 
(ksT <C fi) and the other at higher T. The behaviour 
of the damping rate becomes linear at relatively small 
temperature (fc^T ~ //) in comparison to the homoge- 
neous system B where this regime occurs at ksT 3> \x. 
Moreover, the damping rate increases for larger number 
of condensed atoms because the density of states avail- 
able to the system also increases. It is interesting to note 
that the order of magnitude of the damping rate is the 
same as the one previously estimated for an uniform gas 
and for anisotropic traps [Mpl . 



VI. SUMMARY 

We have considered the monopole oscillation of a Bose- 
condensed dilute atomic gas in an isotropic trap. First 
of all, we have calculated the normal modes of the con- 
densate by solving the time-dependent Gross-Pitaevskii 
equation within Bogoliubov theory |lq and then we have 
used the formalism developed in Ref. [4j to calculate the 
matrix elements associated with the transitions between 
excited states allowed by the monopole selection rules. 
Within a first order perturbation theory we have studied 
the Landau damping of collective modes due to the cou- 
pling with thermal excited levels. We have developed the 
response function formalism to study the fluctuations of 
the system due to an external perturbation. The contri- 
bution of the elementary excitations has been introduced 
also perturbatively as in the calculation of the damping 
strength, and we have derived analytic equations for the 
response function at zero temperature and at low tem- 
perature regime. We have seen that when the condensate 
oscillates with the monopole mode and a small perturba- 
tion to the trap frequency is applied, one can excite new 
resonances at the transition frequencies. These thermal- 
induced resonances are coupled with the monopole due 
to interaction effects. One cannot exclude a priori the 
possibility to observe such resonances also in anisotropic 
traps. This problem deserves further investigation. Ob- 
servation of these resonances would give important and 
unique information about the interaction between ele- 
mentary excitations in Bose-Einstein condensed gases. 
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FIG. 4. Dimensionless damping rate r f/Clu as a function 
of the lorentzian width A (in units of LOho), for No = 50000 
87 Rb atoms in the spherical trap at different temperatures. 
Dashed lines are plotted as a guide for the eye. Solid dots, 
squares and triangles, correspond to the numerical calculation 
at ksT/fi = 0.55, 1.1 and 1.3, respectively. 



FIG. 5. Dimensionless damping rate j/Clu as a function 
of k B T/ii for No = 150000 atoms (dots) and N = 50000 
atoms (squares). Dashed lines are plotted as a guide for the 
eye. 



TABLE I. Damping coefficients 7_r (in units of Who) of 
the "strong resonances" uir (in units of Who) and relative in- 
tensities |j4i/j4o| between temperature-induced the monopole 
resonance, for No — 150000 condensate atoms of 87 Rb in a 
spherical trap with a^a = 0.791 x 10 -4 cm at ksT = /i. 



LUR 



IB. 



\Ai/Ao\ 



1.9115 
2.0252 
2.1432 
2.2576 
2.3655 



0.009268 
0.004147 
0.002102 
0.001097 
0.000545 



0.063 
0.063 
0.127 
0.298 
0.020 



FIG. 1. Transition frequencies w^ (in units of Who) al- 
lowed by the monopole selection rules, for No = 50000 atoms 
of 87 Rb in a spherical trap with a^ — 0.791 x 10~ 4 cm, at 
ksT = /i. The vertical bars have length equal to the cor- 
responding damping strength (in units of wi 10 ). The arrow 
points to the monopole frequency Qm ■ A close-up view of the 
transition frequencies around J7m is presented in the inset in 
order to show the dense background. 



FIG. 2. 

ksT — jj, 



Same as Fig. 1 but for iVo = 5000 atoms at 



FIG. 3. Real part of the response function a(w) (in 
arbitrary units) as a function of w (in units of Who), for 
iV = 150000 at k B T = fi. The dashed line shows the 
monopole resonance at Qm, whereas the other divergences 
of a(w) correspond to the temperature-induced resonances at 
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